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ABSTBiCT 



The specification of the surface heat flux is essential 
for synoptic and seasonal prediction of the upper ocean 
thermal structure. Estimates of -he surface heat flux have 
been prepared for the central Norrh Pacific during January 
1976 through April 1979 using archived fields from rhe Fleet 
Numerical Oceanography Center (FNOC) hemispheric atmospheric 
prediction model. Monthly accumulations of the surface heat 
flux are compared with the change in heat content above 200 
m derived from temperature analyses of the North Pacific Ex- 
periment TRANSPAC ship-of- cpportunit y program. Systematic 
differences are found between the accumulated heat flux 
fields and the oceanic heat content change. Some of the 
differences are due to excessively large changes in ocean 
heat contents above a fixed level. However, our earlier 
studies have suggested a bias of excessive upward surface 
heat flux, especially along the southern boundary of the do- 
main. Assuming local heat balance over a 36-month period, a 
correction field to the FNOC surface heat flux estimates is 
derived. Separate correction fields for the heating and 
cooling seasons demonstrate a seasonal variation in the ac- 
cumulated heat flux versus heat content change values. 
Thus, six bi-monthly correction fields no be added to the 
FNOC heat fluxes are prepared to enable these heat fluxes to 
be used for ocean prediction. 



1, INTR0D0CTI3N 

Prediction of the upper ocean thermal structure re- 
quires a specification of the surface heat fluxes. The time 
scale on which the sensible and latent fluxes, plus the in- 
coming and outgoing radiative fluxes, must be specified is 
dependent on the oceanic phenomena and the type of model. 
Heat flux values resolved on a time scale of 3 hours or less 
must be provided if the diurnal variation in the atmospheric 
forcing is an important consideration in the evolution of 
the upper ocean thermal structure (Garwood, 1977). Consid- 
erable evidence has been accumulated (e.g., Elsberry and 
Camp, 1973; Camp and Elsberry, 197 8; Elsberry and Raney, 
1978; Elsberry and Garwood, 1978; and others) that the upper 
ocean thermal structure responds significantly to atmospher- 
ic synoptic-scale forcing. On seasonal time scales, the 
surface heat flux is accumulated in the upper ocean layers 
during spring and summer, and subsequently removed during 
autumn and winter. 

Ocean mixed layer models provide a means of demonstrat- 
ing the two ways in which the surface heat fluxes affect the 
prediction of ocean thermal structure profiles (Niiler and 
Kraus, 1977). Consider a well-mixed layer of variable 
depth, h, with temperature, T. The heat content (H) per 
unit area in a layer of depth (h) is 

T (2) dz 



where o is the density, Cp is 



(1) 

the specific heat, and z is 
The change in heat content of 



the depth below the surface, 
the mixed layer is 

d (pCpTh) pC),hdT pCpTdh 
dt d t dt 

because the changes in ^ and Cp are relatively small. Solv- 
ing for the temperature change from the first term on the 
right side of (2) gives 

dT 1 d (Th) T dh _ 1 dH 

h dt j)Cj, h dt 



( 2 ) 



dt h dt 

In mixed layer models 
is written as 



(3) 
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entrainment mixing in the Garwood model is dependent on both 
the surface friction velocity and tha surface buoyancy flux, 
which is determined partly by the sic face heat flux. Conse- 
quently, the surface heat flux contributes to both terms in 
(3). The effect in the first term is reflected directly in 
the change in heat content due to the surface flux. The 
thermal structure is strongly dependent on the vertical re- 
distribution of the heat via entrainment mixing, which is 
partly due to the upward surface heat flux. 

We conclude that specification of the surface heat flux 
is an essential factor for prediction of anomalous ocean 
thermal structure. In the next section, we will briefly re- 
view the methods available for estimating the surface heat 
flux over the ocean. In the following section, we examine 
rhe role of the surface heat flux in the oceanic heat budg- 
et. We then derive a correction field to remove a bias in 
the surface heat flux that would be detrimental to ocean 
thermal structure prediction. 



2. SPSCIFICITION OF TH3 S38 FACS H2AT PLOX 
The bulk aerodynamic method for calculating the sensi- 
ble and latent heat fluxes involves the surface wind field 
and the differences in temperature and specific humidity be- 
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tween the sea and the air. Calculations of th= incoming and 
outgoing radiative fluxes require knowledge of the cloud 
cover as well as the solar altitude. An example of a system 
for estimating the surface heat flux from ship observations 
is given by Clark, et al. (1 974), Bunker (1976) and Clark 

(1981). A common approach is to calculare the surface flux 

from each ship report and then average over some space and 
time interval. ?or example, all ship observations taken 
during a month within 5^ latitude and longitude may be used 
to represent the average heat flux in that domain. Husby 
( 1980) has compared the anomalous heat flux over six-month 
periods at Ocean Weather Station V (34*N, 164®S) with the 
estimate based on merchant ship observations within a 4* 
quadrangle. He concludes that the merchant ship reports 
must be carefully screened prior to the heat flux computa- 
tions. 

The accuracy of this method is clearly dependent on the 
number of observations. There is no assurance that the re- 
ports will be randomly distributed in space or time. One 
problem appears to be a "fair weather bias" because the 
ships tend to avoid bad weather. Barnett (1981) indicates 
that a set of heat fluxes in the central Morth Pacific Ocean 
estimated by this method had a bias of 30-45 W m'^ , which 
would result in an excessive estimation of the heat flux 
into the ocean. 

An estimate of the monthly heat flux is not sufficient 
for short-term ocean prediction, which requires information 
on diurnal, or at least synoptic, time scales. Synoptic 
maps of heat flux are difficult to analyze because of the 
sparsity of ship observations. An indirect method is pro- 
posed here. The atmospheric prediction models also require 
a calculation of the heat flux at the ocean surface. The 
heating package for the Fleet Numerical Oceanography Center 
(FNOC) model involves bulk aerodynamic calculations of the 
latent and sensible heat fluxes, plus radiative fluxes that 
are a function of the model-estimated cloudiness (Kesel and 
W inning h off, 19 72) . 
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Our basic hypothesis is that the FNOC atmospheric pre- 
diction model heating package can provide the heat fluxes 
necessary for ocean prediction experiments. Gallacher 
(1979) has described the method used for extracting hourly 
heat flux estimates from the FNOC archives. Elsberry, Gal- 
lacher and Garwood ( 1979) have used these heat fluxes to 
predict the ocean thermal structure changes during the au- 
tumn of 1976. Budd (1980) also used these fields in an at- 
tempt to predict the spring transition from the winter to 
rhe summer regimes in the central North Pacific Ocean. Budd 
found a systematic bias with a too large upward heat flux 
near 30*N. Recently, Steiner (1981) also found a systematic 
bias in these hear flux estimates in the region berween Ha- 
waii and San Francisco. 

The purpose of this paper is to describe a correction 
field to the FNOC surface heat flux fields to permit their 
use in ocean prediction experiments. As in the case of Budd 
(1980) and Steiner (1981) , the approach is to derive a cor- 
rection field that assures a long terra (annual or longer) 
heat balance in the upper ocean. Analyses by Wiite and Bern- 
stein (1979) of the ocean thermal structure observations in 
rhe TRANSPAC ship-of- opportunity program are used to calcu- 
late the time changes in cceanic heat content. These are 
compared with the accumulated surface heat fluxes derived 
from the FNOC files. 5?e then derive the correction to the 
monthly surface heat flux that is necessary to assure local 
heat balance over the 1976-1978 period. In the next section 
we review briefly some tests of the local heat balance as- 
sumption. It is only on the very long time interval that we 
assume local heat balance. We do not require or insure lo- 
cal hear balance on monthly or shorter time intervals over 
which we are doing prediction experiments. 



3. OCEAN HEAT BUDGET STUDIES 
The general purpose in this section is to determine the 
relative importance of the surface heat flux in the upper 
ocean heat budget. Our primary interest is in coen-ccean 
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regimes that would be similar to the conditions in the North 
Pacific Experiment (NORPAX) Anomaly Dynamics Study (ADS) re- 
gion. We do no+ consider regions of boundary currents or 
near-equatorial areas in which horizontal and vertical ad- 
vection are likely to be significant. The space and time 
scales of interest are greater than 1 000 km and one month. 

Many comparisons have been made between the seasonal 
changes of heat content of the upper ocean and the accumu- 
lated heat flux at the surface. Bryan and Schroeder (1960) 
compared the heat content calculated from North Atlantic BT 
data with the surface heating estimated by Budyko (1955). 
They found that the surface heating on a seasonal basis was 
about 20 % less than the change in heat content in the region 
between 20®N and about 50*N. By contrast. Bather. (1971) 
found that the surface heat exchange estimates of Wyrtki 
( 1966) could account for only 29S of the local monthly 
change in heat content in the North Pacific Ocean. Gill and 
Niiler ( 1973) propose that the heat input averaged over 
large areas and times is mainly stored locally, and horizon- 
tal advection by the mean flow is not particularly impor- 
tant. They also cite comparisons by Tabata (1965) and Rob- 
inson (1966) using Ocean Weather Ship data which suggest 
that most of the heat input changes are stored locally. 
Gill and Niiler further suggest that the inaccuracy of the 
heat flux estimates over large areas away from Dcean Weather 
Ships is the likely cause cf some of the departure from lo- 
cal heat balance. 

A careful study of the upper ocean heat budget near OWS 
? based on two weeks cf high-quality observations during the 
Mixed Layer Experiment (MILE) has been reported by Davis, et 
ai. (1981). They found that a one- dimensional upper layer 
heat budget may be closed acceptably if the temperature is 
well-sampled, and if the vertical velocity in the seasonal 
thermocline is also taken into account. 

A recent study by Barnett (19811 used AXBT's along 158** 
and 170*W between 30* and 50*N at ao prcx iraately monthly in- 
tervals to estimate the heat budget. Barnett concludes that 
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approximately 90-95% of the variance in the seasonal change 
of heat storage can be accounted for by the air/sea heat ex- 
change. Horizontal and vertical aivection were of limited 
and no use, respectively, in reproducing the seasonal cycle 
variance. Barnett also concludes that prior heat-budget 
studies that used an arbitrary lower depth in the estimation 
of heat content may have inadvertently included variance on 
different space and time scales. Finally, Barnett suggests 
that the inability to explain adequately the anomalous heat 
content changes in the central Pacific is likely due to im- 
precise knowledge of the source terms in the heat budget. 

Our earlier studies (Elsberry et al. , 1979; Budd, 1980; 
Steiner, 1981) have used an oceanic heat budget to determine 
the feasibility of using the FNOC heat flux estimates for 
ocean prediction. In each case, the imbalances appear to be 
systematic in space and time. Large, short-term imbalances 
are likely due to ocean observational sampling errors and to 
using an arbitrary lower depth in the estimate of heat con- 
tent (Barnett, 1981). However, the long-term, systematic 
error is likely due to a bias in the FHOC heat flux esti- 
mates. In the following sections, we derive an appropriate 
correction field to be added to the FNOC heat flux estimates 
to remove the long term bias. 

4. CALCDIATION CF THE OCEAN HEAT BUDGET 

Monthly mean temperatures during 1976-1979 have been 
objectively analyzed by White and Bernstein ( 1979) based on 
TRANSPAC shi?-of-opportun i ty XBT’s only. Analyzed tempera- 
tures are available on a 2* latitude by 5 ^ longitude grid at 
0, 20, 40, 60, SO, 120, 1 50, 200, 250 , 300 and 400 m. The 

general domain of the heat budget calculations is from 30* 
to 50*N and from 170*E to 130®W. However, the number of ob- 
servations in the northwest and southwest corners of the 
grid are insufficient to provide reliable estimates during 
some months. Consequently, these areas are eliminated from 
the following analysis, and will appear as cross-hatched 
areas on all figures. 
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To calculate the heat content as in (1), the objective- 
ly analyzed temperatures are first interpolated to 5 m 
depths. The heat content is actually calculated relative to 
the 200 m temperature, with the intent of renoving some of 
the change in heat content due to vertical die placements of 
the thermccline. Barnett (1981) has suggested that selec- 
tion of an arbitrary depth such as 20 0 m may introduce addi- 
tional variance into the heat content calculation due to in- 
ternal waves, tides and so forth. Inspection of 

month-to-month heat content values does reveal time changes 
which are much larger than can be reasonably accounted for 
by surface heat fluxes. Large positive time changes tend to 
be followed by large negative changes. Thus, time intervals 
of either two or three months have been tested in the esti- 
mate of the oceanic heat content change, 

= H(t + ^t) - H (t) . (5) 

We expect to be positive during the period of net down- 

ward surface heat flux (roughly between April and September) 
and negative during the period of net upward flux. The 
monthly T5ANSPAC analyses are assumed to apply on the 15th 
of each month. If ^t is two months in (5) , then will 

also apply at the 15th of the month. For example, the dif- 
ference between the March 1976 and the January 1976 heat 
content fields would be applied at the middle of February. 
Notice that this is the first possible value that can be 
calculated because the analyses begin in January 1976. How- 
ever, if ^t is three months in (5), then AH will apply at 
the mid-pcint of this interval. Thus, the first possible 
difference would be between April 1976 and January 1976, and 
it would apply on 1 March 1976. The advantage of the longer 
time interval is that shorter period fluctuations in heat 
content tend to be averaged out and one obtains a more con- 
sistent measure of the seasonal variation. The disadvantage 
of the longer time interval is the greater inaccuracy in the 
finite difference approximation to i derivative. Finally, 
it is necessary to extend the AH calculation into 1979 to 
obtain three complete annual cycles. In the case of A" ' 2 
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(3) months, the analyses through February (March) 1979 are 
used. 



Following Emery (1976), 
may be written as 

^ - ^(H - jCpDT ) 
^ ~ D 



the basic heat balance equation 



- V • + Q , 



( 6 ) 



where H is defined as in (1) with -h , the mixed layer 
depth, replaced by D which is a (constant) depth below the 
main thermocline. In (6) , the first term or. the right rep- 
resents the divergence of heat due to horizontal divergence, 
WpH/D, and vertical advection The second term is 

the horizontal advection cf heat and the last term is rhe 
surface heat flux. The suggestion by Barnett (1981) and 
Davis et . al. (1981) that D be the depth of an isorherm be- 
low the main thermocline rather than a constant depth re- 
sults in the divergence term in (6) being zero, and thus 
eliminates a source of uncertainty from the heat budget. 
Due to the large vertical sampling interval in the TRANSPAC 
analyses, we will use D = 200m throughout this study. Thus, 
we assume the following budget equation 

= j Qdt Residual = Q-. + Residual . (7) 

The integral sign indicates that the air-sea flux is summed 
over the same time interval that ^8 is evaluated. Q.j. is 
defined as che surface heating. In the ^t = 2 month exam- 
ple above, the sum would be from 15 January through 15 
March. The residual term in (7) includes the remaining 
terms (physical effects) in (6) that can not be evaluated, 
plus the errors in estimating the heat content changes and 
the surface fluxes. If these observational errors are ran- 
dom over a sufficiently long interval, their contribution 
should be averaged cut. In particular, we sum over 36 
months (or a small subset if a seasonally varying correction 
factor is derived) and derive an averaged correction field 
to be added to the surface heat flux. However, the physical 
terms may not sum to zero. Persistent vertical and horizon- 
tal advection or diffusive effects that accumulate during 
the period are thus also incorporated into the correction 
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field. The impact of assuming a local (one-dimensional) 
heat balance via this long-term correction field will be as- 
sessed in the prediction experiments. 

The total heat flux Q in (7) may be expressed as 

C = - (Qb * * ^1^ ' 

where the subscripts s, b, 1 and h refer to solar, back, la- 
tent and sensible heat fluxes through the sea surface. A 
description of the calculation of each of these heat fluxes 
in the FNOC atmospheric prediction model is given in Kesel 
and Minninghoff (1972). The values of these terms after a 
one-hour integration of the model are taken to apply at the 
synoptic map time (00 or 12 GMT). After spatial interpola- 
tion to the grid points of the ocean analysis, these instan- 
taneous values are interpolated to hourly values using tech- 
niques described by Gallacher (1979). The representation of 
the solar flux is particularly involved. It is this time 
series of hourly values of and Qj, (= that 

are accumulated in (7). Simply summing the instantaneous 
12-h FNOC values would not produce the same sum because we 
have reconstructed the hourly variation in Q. 

An example of a field of near the beginning of the 

ocean heating season is given in Fig. 1 . Over most (all) 
of the domain, the heat content change is positive, as ex- 
pected. A value of 2 x 10 cal cm over this two-month 

period corresponds to an average inr rease in temperature of 
2® C over a depth of 100 m, or over a depth of 50 m. 

Less confidence should be given to the values near the 
northern and southern boundaries of the domain. The primary 
ship-cf- opportunity tracks lie on the great rircle route be- 
tween Japan and the USA west coast ports ('rfhite and Bern- 
stein, 1 979). During the winter, the tracks tend to be dis- 
placed southward, which makes the northern region sampling 
rather poor. Somewhat the opposite effect occurs during the 
summer as the ship tracks are extended northward. The ob- 
iective analysis procedure (Shite and Bernstein, 1979) used 
thrcuah 1 978 would generate values at all gridpoints regard- 
less of whether any observations fell in the vicinity of the 
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point. Because the deficient points tend to lie toward the 
north and the south, the fictitious values tend to resemble 
northward or southward extrapolations from the central re- 
gion. 

The corresponding integrated surface heat flux for this 
period derived from the FNOC calculations is shown in Fig. 
1c. In contrast to the fields, no extrapolation of the 

fields toward the domain boundaries is involved in the 
field. In the band from UO-50® N, the surface exchange is 
approximately equal to the observed ocean heat content 
change, which suggests an approximate local heat balance. 
However, the net FNOC heat flux continues to be upward be- 
tween 30-35*® N during this period. One can not tell from 
this diagram whether the daily values of in (8) are too 
small, or if the upward surface heat flux (Q^,) is too large. 
The fact that the maximum values are found near the longi- 
tude for which the boundary of the FNOC grid is tangent to 
the equator (therefore closest to the study region) suggests 
that the problem may lie in the boundary conditions that are 
a pplied . 

The difference between and the cumulative heat flux 

(Q^) is shown in Fig. 1 . A light filter, with weights 
2-4-2 in latitude and 1-4-1 in longitude, has been passed 
over this field to reduce small scale noise. Although there 
are considerable areas with the expected zero values, there 
are other areas with positive and negative values. In par- 
ticular, the residual term in (8) is very large along the 
southern boundary because of the lack of a net downward sur- 
face heat flux (Fig. 1 ) . 

The complete set of monthly , 2^ and values 

is contained in the Appendix. Since our desire is to obtain 
a correction field, it is not appropriate to discuss hare 
each of these sets. One general feature is that the 
fields have less east-west variation and more north-south 
variation than do the fields. There are several possi- 
ble causes for this feature: (a) the FNOC heat fluxes may 

not contain an adequate renresentation of the east-west var- 
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iations in the anomalous surface fluxes; (b) vertical advec- 
tive processes associated with ocean eddies are contained in 
the heat content values relative to a fixed depth; and (c) 
errors due to inadequate observations or analysis of the 
ocean heat content. The spring transition regime example in 
Fig. 1 illustrates one of the worst imbalances in the time 
series. An excessive upward heat flux also occurs along the 
southern boundary during autumn. However, the H field is 
negative (upper ocean heat content diminishing in time) dur- 
ing this period. Consequently, the magnitude and sign of 
the imbalance may have a seasonal dependence (discussed fur- 
ther in the next section). 

Similar sets of A H, Q-,. and A H-Q^ were also prepared 
with ^ months. These fields (not shown) are not di- 

rectly comparable with the two-month set because of nhe 
half-month displacement of the central point. However, the 
time evolution of the major features using A- == 3 months 

can be easily associated with the evolution shown by the ^t 
= 2 month series. Therefore, only the correction fields 
based on the 2-month differences will be discussed below. 

5. INTZRPSETATION OP THE HEAT BUDGET IMBALANCE 

Perhaps the most important factor to emphasize regard- 
ing these calculations of the local heat budget (8) is that 
each of the two terms is subject to large errors. As noced 
above, calculating the heat content changes relative to a 
fixed depth may produce oscillations equivalent to a 1-2^C 
change averaged over the entire depth. A faulty XBT could 
also be the cause of a temperature bias with depth. How- 
ever, we would expect that the objective analysis technique 
would tend to eliminate such a bias if other correct pro- 
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Because the heat budget imbalances may arise from ei- 
ther term. Table 1 was prepared to indicate tha various pos- 
sibilities that may arise. It is useful to separate the AH 
into separate periods when the ocean heat content is in- 
creasing cr decreasing. Apart from small calculation er- 
rors, one would hope that the corresponding would be po- 
sitive (surface heating) or negative, respectively, so that 
Cases A-B or E-F would apply. One expects approximately an 
equal distribution between A and B or between E and F if the 

errors are random. The percentage of gridpoints with month- 

11 -2 

ly <AH~Qx differences exceeding 0.5 x 10 cal cm are nor 
evenly distributed between cases A and 3 or between E and F 
when the uncorrecred is used in the differences. There 



is clearly a bias toward cases A and 



which could be at- 



tributed to excessive upward heat flux during both the ocean 
warming and cooling periods. However an approximately equal 
distribution is obtained when the corrected (see description 
of six bi-monthly correction fields in the next section) 
surface heat flux is used in the differences. 

A value of 0.5 x 10 cal cm corresponds to a heat flux 
bias of 7.0 cal cm^hr*or a temperature bias of ,25* C in a 
200 m water column. Both of these values are within the ex- 
pected range of instrument error for the measurements which 
were used in the analyses. Whereas 6535 of all the differ- 



ences exceed this criterion for the uncorrected heat 
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only 42'’o exceed the criterion for the corrected heat flux. 

Cases C and D in Table 1 are labelled as drastic imba- 
lances because the and are of opposite tendencies. A 

physical explanation for such an event might be the advec- 
tion of a warm (cold) ocean eddy into the region that has 
upward (downward) surface heat flux. An example of a compu- 
tational explanation is found along the southern boundary in 
Fig. 1 where < 0 when > 0. The percentages in Table 

1 indicate that drastic imbalances are relatiyely rare dur- 
ing ocean cooling periods (Case D - 1%) when the uncorrected 
is used in the differences. However, this is not the 



rase during the ocean warming periods. 
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wh en 



fluxes rathec 



of the points have significant upward heat 
than downward. Since the heat budget for the central Pacif- 
ic is in an approximate local balance for seasonal time 
scales (Gill and Niiler, 1973), the percentages for Cases C 
and D should be less than those for Cases A, B, E and F, 
This distribution is obtained when the corrected values 

are used, but nor with the uncorrected 3^. This is a fur- 
ther demonstration of the need to adjust the FNOC surface 
heat fluxes. 

The alternative of a systematic bias due to the AH 
calculation does not appear likely. Although there are more 
points with AH positive (59%) than negative this as- 

ymmetry is probably within the limits of roughly offsetting 
periods of ocean warming versus cooling. The imbalances be- 
tween seasons in Table 1 are not consistent with horizontal 
advection being a primary cause. One would expect stronger 
Ekman advection effects during the winter (cooling season), 
whereas the larger imbalances are found during the summer. 



TABLE 1 

A H versus 2^ 





^ ^ . 

(Downward surraca 
heat flux) 


Qt < 0 

(Upward surface 
"heat flux) 


AH > 0. 
Ocean warmingi 


Oase A /\ H — Q y ^0 

Inadequate downward Q 
Excessive ocean warming 
22.4% (26.8%) 


Case C A^H— 2x ^ ^ 

Drastic imbaian.ee 
33.4% (7.3%) 


Case B AH-Qt < 0 

Excessive downward Q 
Inadequate ocean warming 
2.9% (20.0%) 


AH < 0 
Ocean coding 


Case D AH-Qt < 0 

Drastic imbalance 
1.0% (6.5%) 


Case E AiH-Qr > 0 

Excessive uoward Q 
Inadequate ocean coolin 
34.g% (16.7%) 


Case ? AH-Qt < 0 

Inadequate upward 0 
Excessive ocean cooling 
5.4% (22.6%) 



Possible ohysicai association for oositive or neaatrve neat buacet 
imbalances aiven that the heat content or integrated surface rluxes 
are positive cr necative. Percentage of the gfidooints for the 36 
monthly maos that fall in each category are ihiicated. The o<^rcen*aae 
in parentheses were ccmouted using th e' corrected total heat' flux. 

Only points with 1A'H“QtI ^ 0-5 x 10'*' cal cm'* for the two month 
difrerence are inc-udea in the calculation of the oercent. 



6. CALCOLATIOH OF CORRSCTIOM FIELD 

Our objective is to determine a correction field that 
may be applied to the Q values in (9) during the ocean pre- 
diction experiments. Because the time step in the Garwood 
model is typically 1 h, it is desirable that the correction 
be in units of cal cm'^ hour”*. This is done by converting 
each of fields to an hourly rate using the proper 

number of hours in the time interval. The sign is also 
changed because the correction is to be added to the Qy 
f i eld , 

It is not within the scope of this paper to determine 
the fraction of this correction which will be applied to Qj 
versus the remaining -hree terms in (8). It might be noted 
that occurs only during the daytime and this flux is to 
be distributed exponentially with depth by the prediction 
model. In contrast, the remaining three terms in (8) apply 
only at the surface and tend to be positive (upward) 
throughout the day. Thus, one may expect considerable dif- 
ferences in the upper ocean predictions as the fraction of 
the correction that is allcted to is increased. 

For simplicity, one desires a single correction field 
as in Fig. 2. The basic features in this correction field 
are generally consistent with the pattern shown in Fig. la. 
In particular, the correction reduces the surface heat flux 
(upward is defined to be positive) along the southern bound- 
ary. It can be shown that this correction is equivalent to 
calculating the residual in (7) using the difference in H 
between Januar y-February 1979 and Jan uary- February 1976, and 
the Qy for the entire period. He may also regard the single 
correction field as the adjustment necessary during each 
hourly time step of an integration from 15 January 1976 to 
15 January 1979. If the correction is applied in this way, 
we can be assured of conservation of heat at each gridpoint 
for the entire three year period (assuming no vertical dif- 
fusion and no round-off errors) . 

Each of the 36 maps of is different from that 

implied by the single correction field. He noted above that 
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the seasonal variation in or night be associated with 

a modulation in the difference field. A separate correction 
field is derived for November through April (Fig. 3a) and 
for May through October (Fig. 3b). Each of these correc- 
tions is summed over the three years and six months, so tha- 
a total of 18 values is included in each average. There is 
clearly a seasonal variation in these two correction fields. 
The pattern during the heating season (Fig. 3b) is similar 
to the single correction field (Fig. 2), but the values are 
larger. The pattern during the cooling season (Fig. 3a) is 
quite different. Not only are the values considerably 
smaller than in Fig. 3b, there is more of a north-south ori- 
entation of the isolines in the north-central region. 

Based on the seasonality in Fig. 3, we also examined 
further subdivisions into quarterly or bi-monthly correction 
fields. A decision was then made to adopt the six correc- 
tion fields in Fig. 4 a-f. It is felt that six maps will 
give a better representat icn of the seasonal variation. In- 
spection of Fig. 4 b-e reveals a basically east-west pat- 
tern, whereas the remaining two maps exhibit the north-south 
orientation. Even though the seasonal variation does not 
appear to be a sine wave, the six maps appear to provide a 
relatively smooth transition between the two basic patterns 
in Fig. 3. There are some non-sea sonal features in these 
correction fields. The largest is found along the western 
side of the domain. Larger corrections are required during 
February- Marc h (Fig. 4a) and during August-September (Fig. 
4d). The transition between the correction field for Decera- 
ber-January (Fig. 4e) and for Feb- a ary-March is especially 
noteworthy along the western boundary. The regions of posi- 
tive corrections (additional upward heat flux required) also 
tend to be somewhat erratic. However, these positive cor- 
rections are always small, so they will have little effect 
on the ocean predictions. 

An example of the application of the heat flux correc- 
tion field at a specific location is given in Fig. 5. The 
seasonal changes in neat content calculated from th- 
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TRANSP&C analyses have rslatively constant values during the 
heating season. By contrast, the values during the 

cooling season tend to reach rather well-defined peaks. The 
surface heating agrees rather well with the maximum val- 

ues during each coding season. However, the uncorrected 
heat flux during the warming season is clearly deficient, 
especially during 1978. Applying the heat flax corrections 
for this location from Fig. 4 improves the agreement between 
the surface heating and the H. The major feature to be 
noticed is that the H - Q^) differences are indeed sys- 

remaric, so that a single correction in each two-month peri- 
od tends to improve all three years. One finds periods in 
which the - Q.p remains large after the correction has 

been applied. This indicates that although we have had to 
correct the surface heat flux we have not forced the heat 
budget to be one-dimensional for periods shorter than 36 
months. Some examples are the differences during April - 
June 1976 and during July - August 1978. Another feature in 
Fig, 5 is that the area between the and the corrected 

curves must sum to the difference between the final and ini- 
tial H values (February 1 979 and January 1976, respective- 
ly) . This requirement is a consequence of the local heat 
budget assumption over the 36 -month period. 

7, SOMMAHI 

We have prepared six bi-monthly correction fields to be 
applied to the FNOC heat flux values to be used for ocean 
prediction. The largest corrections are found generally be- 
tween 30 and 38 N during April through November. The heat 
flux bias is evidently not serious for the acmospheric pre- 
dictions because they are limited to 72 h . However, such a 
bias can be disastrous (Budd, 1980> for ocsan prediction 
over monthly time periods. The correction fields are aver- 
aged over three years (1976-1 979) . Closure of the local 
heat budget over the entire period is insured by the use of 
this correction field. However, this is not true over any 
shorter time intervals. In particular, fluctuations in heat 
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content in the monthly TRANSPAC analyses are not accounted 
for by the long-term correction field. 

The success of these correction fields can only be 
judged by their application in the ocean prediction experi- 
ments. Inclusion of the corrections should result in im- 
proved predictions. One of the tests will be whether the 
six bi-monthly correction fields perform better than the two 
semi-annual correction fields. In some of the correction 
fields, the east-west variation is not large. It is possi- 
ble that a correction dependent only on latitude may perform 
as well. Such a correction would be easier to apply. The 
numerical ocean prediction experiments necessary to demon- 
strate the usefulness of the corrections are in progress, 
and will be reported separately. 
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LIST OF FICaRBS 



Fig. 1 The difference (A change in heat 

content relative to 200 m (AH) and the surface heat- 
ing Q, ) between July 1976 and May 1976. Cross- 
haxched areas have insufficient data for xhe AH 
analysis. Positive values of AH indicate increasing 
heat content, while posixive Q- indicates a net down- 
ward surface heat flux. Units: 10' cal cm** oer 2 

months. 

-i - j 

Fio. 2 Correction field (cal cm h | to be applied 
to FNOC heat flux fields based on all 36 of the 
monthly sens of j(AH-0») evaluaxed over two-month in- 
tervals. Negative values indicate that the upward 
heat flux is to be reduced by the amount shown. 

Fig- 3 As in Fig. 2, except separate correction 
fields for November through April (topi and for May 
through October (bottom). 

Fig. 4 As in Fig. 2, except bi-monthly correction 
fields for: February - March 'Ctop) ; April - May (mid- 
dle) ; June - July (tottom). 

Fig. 4 (continued) As in Fig. 2, excaot bi-monthly 
corrections for: August - September (top); October - 

November (middle) ; December - January (bottom) . 

Fiq. 5 Time series of changes in heat content 
(dashed), uncorrected surface heating (solid)_, and 
corrected surface heating (dotted) ax 35®N, 175®W. CQ, 
is the corrected total surrace heating. Units are 10' 
cal cm"^ over two- month intervals. 



SON 




Fig. 1 The difference change in heat 

content relative to 200 ai (A-^) ^nd the surface heat- 
ing (Q, ) between July 19 /b and ilay 1976. Cross- 
hatche'd areas have insufficient data for the AH 
analysis. Positive values of indrcate increasing 

heat content, while oositive indicates a net down- 
ward surface heat flux. Units: 10' cal cm*^ per 2 

nonths. 
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Fiq. 2 Correction field (cal cm h ) to be applied 
to FNOC heat flux fields based on all 36 of the 
monthly sets of j(A H-Q,- 1 evaluated over two-month in- 
tervals. Negative values indicate that the upward 
heat flux is to be reduced by the amount shown. 
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Pig. 3 As in Fig. 2, excapt separate correction 
fields for November through April (tppt and for May 
through October (bottom). 
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4 As in Fig, 2 , except bi-monthlv correct 
for: February - March (top) ; April - May (m 
June - July (bortom) . 
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Fiq. 5 Time series of ohangs 
(dashed) ^ uncorrectea surface ne 
correctea surface heating (dotted) 
is the corrected total surrace hea 
cal cm"^ over two-mcnth intervals. 
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Appendix A 

MONTHLY SETS OF H AND 

The monthly sets of ^ H and Qy based on two- 

monthly intervals are included, except for the Hay - July 
1 976 set which is in Fig. 1. The form of these diagrams is 
described in the caption of Fig. 1. Positive values (solid 
lines) are increasing heat content in time and net surface 
heating, and conversely fcr negative values (dashed lines). 
Zero lines are enhanced. 
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